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^ ■ Abstract 

o : 

O . Cobalt nanowires with a diameter in the range between 50 to 

lOOnm can be prepared as single-crystal wires with the easy axis 
(the c-axis) perpendicular to the wire axis. The competition between 
, the crystal anisotropy and demagnetization energy frustrates the mag- 

\^ I netization direction. A periodic modulation of the angle 6 between M 

' and the wire axis yields a lower energy. 

P : PACS: 

^ ■ In recent years magnetic nanostructures have eperienced a great interest 

O ■ in dynamic magnetic torque experiments (for ref. see for example [1]). These 

experiments explore the possibility to rotate the orientation of the magneti- 
^ [ zation with a current pulse. This could be an important tool in spintronics. 

5^ I Complementary to the dynamic experiments we want to explore the static 

properties of magnetic nanostructures, in particular cobalt nanowires (NW). 
We believe that a detailed knowledge of the static magnetic properties will 
have important consequences in their dynamic behavior. 

A number of experimental groups f2], [3], [4], [5], [6], [7], [8], [9] have pre- 
pared Co NWs with diameters in the range of 30nm to several lOOnm. Similar 
Co NWs with a diameter of 80nm were recently fabricated at the University 
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of Southern California [TO]. In some of the experiments the magnetic struc- 
ture of the NWs was investigated with a magnetic force microscope (MFM) 
[2], [8], |7], [To]. The MFM scan showed spacial oscillations of the magnetic 
field along the length of the wire which are sometimes quasi-periodic. Thiav- 
ille et al. [11] concluded that in their experiments the period is in agreement 
with a "head-to-head" magnetization, a 180° Bloch wall. 

Henry at al. [7j observed by means of dark field transmission electron 
microscopy (TEM) images that the Co NWs have the bulk hexagonal struc- 
ture. For wire diameters 2R < 50nm the easy c-axis lies parallel to the wire 
axis while for NWs with diameters of 2R > 50nm the easy c-axis is perpen- 
dicular to the NW axis. In the following we discuss the latter case, NWs 
with a diameter 2R > 50nm. Below we choose a radius of R = 40nm for 
quantitative calculations. We denote the wire axis as the z-direction and the 
easy axis as the x-direction of our coordinate system. 

When the axis of the Co NW and the easy axis lie perpendicular to each 
other then the magnetization is frustrated. The shape or demagnetization 
anisotropy prefers to align the magnetization in the z-direction, parallel to 
the wire axis. But the crystal anisotropy definitely favors the x-direction. 
And this crystal anisotropy is very large in the uniaxial Co wire. 

The shape or demagnetization anisotropy energy density (ED) is due to 
the demagnetization field and given by 

2 2 

where N is the 3x3 demagnetization matrix, H is the demagnetization field, 
i.e. the magnetic field in the absence of an external magnetic field, and M 
the magnetization. We introduce 6 as the angle between the z- or wire axis 
and the magnetization M. (Within this paper the magnetization will always 
lie in the x-z-plane). Then one has a demagnetization factor of N^x = \ for 
^ = 7r/2 (perpendicular to the wire) and N^z = for 6* = (parallel to the 
wire axis). For a constant magnetization M = Mo (sin 6*, 0, cos 6*) under the 
angle 9 the demagnetization energy density is 



We take from O'Handley [12] the reference value for the magnetic ED of Co 
= ifM^ = 12 X lOV/m^. The value of Wd/uoo is for M parallel to the 
wire axis and 1/2 for M parallel to the easy axis. 
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The energy density of the crystal anisotropy is generally given in terms 
of the angle between the easy axis and the magnetization. In our geometry 
this angle is equal to (f — 6') . The crystal anisotropy ED is, in terms of this 
angle 6 

Uca = ki cos^ 9 + k2 cos"^ 9 

The crystal anisotropy constant ki is given in the literature as fci = 4.1 x 
10^ J /m^ [12]. For the constant k2 one finds different values in the literature, 
for example fca = 1.5 x lOV/m^ [12j and ^2 = 1.0 x lOV/m^ [13]. The 
resulting crystal anisotropy ED Uca/uoo is 0.47 (0.425) along the wire axis 
and parallel to the easy axis. (The value in parenthesis is for k2 = 1.0 x 
10^ J/m^). This difference in the constant /c2 has important consequences. 
The value of ^2 = 1.5 x 10^ J/m^ yields the lowest ED {ud + Uca) /uoo = 0.450 
for a finite angle of = 0.65 = 37° between the magnetization and the wire 
axis. For the value of k2 = 1.0 x 10^ J/m^ the magnetization would align 
parallel to the z-axis. 

Obviously the competition between the crystal anisotropy and demagne- 
tization is a close call. The system will try to reduce its energy as much as 
possible by the crystal anisotropy without paying too much energy to the 
demagnetization energy. One way to reduce the demagnetization energy is 
to modulate the magnetization direction in the x-z-plane so that the angle 
9 between M and z oscillates as ^ = ^qCos {kz). (There is no oscillation in 
time but only in space in contrast to spin waves in NWs which have been 
treated by Arias and Mills p^). While for a constant magnetization in x- 
direction the field H falls off as 1/ with the distance p from the wire axis, 
a modulated magnetization with a period A will cancel the field for distances 
p which are larger than A. This reduces the demagnetization ED. In this 
paper we investigate the effect of such a modulation on the ED of the wire. 
This modulation corresponds to a magnetization M 

M = Mq (sin (6*0 cos kz) , 0, cos (6*0 cos kz)) (1) 

In Fig.l the orientation of the magnetization is shown as a function of z. We 
keep the absolute value of |M| = Mq constant. 
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Fig.l: Nanowire with magnetization modulation along the z-axis. 



The magnetization components and can be expressed as two 
Fourier series. 

{z) = MoE"=oC2.+i cos [{2u + 1) kz] 
M, {z) = MoEr=iC2. cos {2ukz) 

The coefficients C2u+i, C2u can be easily obtained from a Fourier expansion of 
M in equ. ([1]) . The lowest coefficients are Cq (6*0) = (l — 16*0 + -^0^ — h..) , 
ci (^0) = (^0 — 1^0 + '^^0 — '"••)' ^^c. We include terms up to the order of 

In the next step we calculate the demagnetization field H for a magne- 
tization Mj; = Mj;o cos (g^;) . Setting afterwards q = (2z/ + 1) k and M^o = 
M0C2U+1 the results can be used for each Fourier component. 

The magnetic flux B inside and outside of the sample is given by B = 
/io (H + M). Since there are no external currents in our problem the curl of 
the magnetic field vanishes, V x H = 0. Therefore the magnetic field can be 
expressed as the gradient of a magnetic potential H = — V0 (in full analogy 
the electrostatic case). Taking the divergence of the magnetic flux (which 
vanishes) yields 

= V ■ B = /io(V ■ H+V ■ M) 
and replacing the field by the potential yields 

A0=V ■ M 

For Mx the divergence of M is zero. 

We use cylindrical coordinates (p, ip, z) and take the (/^-dependence as 
cos if. Then the solutions of the LaPlace equation are 

, _ J (gp) cos if cos qz p < i? 1 

I C°'^^Ki (qp) cos (f cos qz p > R j 
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where Ji (s) and Ki (s) are modified Bessel functions. The coefficients C", C°"* 
are obtained by using the boundary conditions ai p = R. The components 
Bp and if^ have to be continuous. This yields C*" = RM^qKi (qR) and 
(jout _ jiji^^^j^ (qR)- (In determing the coefficients one obtains the Wron- 
ski determinant W = [h (qR) K[ (qR) — I[ (qR) Ki (qR)] as a denominator, 
which has the value W — —1/ (qR))- 

From the magnetic potential one obtains the components of the magnetic 
field H. The x-component of H inside the wire is 



{p<R) = -qRK, {qR) 



I[ (qp) cos^ (fi + —Ii(qp) sin^ ip 



Mq^ cos (qz) 



The local demagnetization ED is — {po/2) H^M^. We average over a period 
in z-direction and the cross section ttR"^ and obtain for an individual Fourier 
component the demagnetization ED 



( yM^^o) (qR) h {qR) 



For each q = (2z/ + 1) k the demagnetization field H interacts only with 
the magnetization M of the same q (after averaging). Then the total contri- 
bution of all components of is just the sum of the individual contributions. 
In the following we normalize all EDs by dividing by the value uqq = ^Mq. 
Then the normalized ED is 

= ELo \k, [{2v + 1) .] h + 1) .] 

Moo ^ 

where s = kR. In the numerical evaluation we include three terms (the third 
hardly contributes). 

The Fourier components — M^o cos qz for the z-component of the 
magnetization are calculated quite analogously. The main difference is that 
the magnetic field H and therefore the magnetic potential are independent of 
(f. Therefore (f) is given by the modified Bessel functions Iq (qp) and Kq (qp) 
of order zero. Furthermore A0 does not vanish but is given by 

A(/'=^ = -M^ogsm (qz) ^ 

The solution is found in complete analogy to the M^-component and is 
given by 
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(j) (p, z) = RMzo sin (qz) 



^ + K',{qR)Io{qp)\ p<R 
I'^ {qR) Ko (qp) p>R 



The magnetic field component inside the wire is 

H,{p<R) = - {qRK'o {qR) Jq (qp) + 1) Mq,, cos qz 

In the evaluation of the demagnetization ED we use the identities Hq (t) = 
dithit)) /dt, K'^it) = -Ki{t), r^^it) = hit). The averaged demagneti- 
zation ED becomes (^Mj^) (| - {qR) {qR)). The contribution of all 
Fourier components of is 



Uz {s,9o) 



Moo 



Er=i (C2. {e))' Q - K, {2us) h (2^.)) 



Again we include the first three terms in the numerical evaluation. 

Next we consider the crystal anisotropy ED. The average of the term 
ki cos^ 6 yields 

(1) (aA 11 /-stt 



u)^ {9o) 



1 1 f 

= / ki cos^ (610 cos {s)) d {s) = 0.34 x ai {Oq) 



where ai(^o) = l~|^o + |^o — ^he average of the term k2 cos^ 9 yields 



(2) 

8. 3 X 10"^ X as (^o) 



„ fa \ — 1 a2 , 



for k2 ^ 1.5 X 10^ J/m^ with 02 (^0) = 1 - ^0 + f^d ^ +■■■ both cases we 
include terms up to the order of ^^q^. 

Finally we have to include the exchange stiffness of the Co wire. While 
a modulation of the magnetization can reduce the demagnetization and the 
crystal anisotropy energy, it will cost energy because the of the bending of 
the magnetization. The increase in the ED can be expressed in terms of the 
exchange stiffness constant D^x 

u --^e'D k' 
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Since the energy densities and function of kR and 6q we express 

all energies as functions of s = kR and 9q. Then we obtain 

Uex _ 1 MqD^x q2 ^2 
Moo 'iUooR'^gfiB ° 

Liu et al. [15] determined the exchange stiffness D^x experimentally from 
the spin-wave spectrum in hexagonal Co. They also performed a theoretical 
calculation. From the experiment they obtained Dg^ = 4:35meV x = 
6.96 X 10~^° Jm^. Their theoretical result yielded twice this value. Using 
the experimental value and a radius of R = 40nm we obtain u^x = 8125 
* [kR)'^ 6q [Jm~^]. The normalized exchange stiffness ED is then 

— = a^xS^el, a^x = 0.6 8 X 10"^ 
Moo 

This exchange ED is very small compared with the demagnetization and the 
crystal anisotropy EDs which are of the order of 1. 

Finally we add all terms and calculate the total ED as a function of 
s = kR and 6o and determine the minimum of this energy 

Ut (S, 0) = — [Ux {S, 9o) + (S, Oq) + Uca (Oq) + U^x (s, Oq)] 

We perform the calculation for different choices of the parameter k2 and 
determine the position of the minimum of the ED in the s-Oo-plane. To 
investigate the effect of the exchange ED we also perform a calculation with 
twice the experimental value for aex- In table I the numerical results for 
different parameter choices are collected. 



k2/M00 


Biex/Uoo [10-^] 


Smin 


^min 


Umin/ Uqo 


u[M||z] 


u[M||X] 





0.6 8 


2.3 


0.7 


0.333 88 


0.342 


.5 





1.36 


1.75 


0.3 


0.34137 


0.342 


.5 


.083 


0.6 8 


2.1 


1.0 


0.378 83 


0.425 


.5 


.083 


1.36 


1.6 


0.8 


0.396 9 


0.425 


.5 


0.125 


0.68 


2.1 


1.0 


0.39704 


0.425 


.5 


0.125 


1.36 


1.5 


0.9 


0.417 71 


0.425 


.5 



Table I: For two values of ^2 and a^x the coordinates and the value of the 
(normalized) energy density (ED) in the s — ^^o-plane are collected in columns 
three, four and five. Columns six and seven give the ED for a constant 
magnization parallel to the z- and the x-asis. 
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For k^/uQQ = 0.125 and a^x = 0.68 x 10~^ we find the minimum at 
{s,6q) = (2.1, 1.0). In Fig.2a,b the dependence of Uf/uQQ is plotted for these 
parameters. The figures show two orthogonal traces through the energy 
minimum (a) along the s = kR direction and (b) along the direction. 



0.55 T 




Fig.2a,b: The ED as a function of s = kR (3a) and Oq (3b) through 

the minimum for the parameters k'z/uQQ — 0.125 and Oex/i^oo — 0.6 8 x 10 

For ^0 = 1-0 we can draw the two components and as a function 
of z along the wire. This is shown in Fig. 3. With 9 — 1.0 * cos(A;2;) the 
amplitude of the angle is less than 7r/2. Therefore the z-component never 
reverses direction. At sin (1.1) = 0.84 the x-component reaches almost the 
saturation magnetization. 



1.0-, 




Fig. 3a: The x- and the z-component of the 
magnetization as a function of position s — kz. 
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For a comparison we calcTilatc the ED when the magnetization angle 
rotates in the x-z-plane as 6' = cos (kz). This yields 

M =Mo (sin (kz) , 0, cos (kz)) 

In this case we have only one Fourier component in x- and z-directions with 
the same wave number k. The demagnetization ED follows from the above 
calculation. (There is no cross term between the x- and z-part of the de- 
magnetization ED since their 99-components are orthogonal). The fci-part 
of the crystal anisotropy ED has the weight 1/2 and the k2 has a weight of 
3/8. The exchange stiffness ED is just Uex/uoo = 0.68 x lO^^s^. Fig.4 shows 
the dependence of the total ED Ut/uoo as a function of s = kR. The total 
ED has its smallest value of Ut/uoo — 0.4677 at A; = 0. This value is con- 
siderably higher than for a constant magnetization along the NW axis with 
u\\/uoo = 0.408. Therefore this behavior of the magnetization is energetically 
unfavorable. 

0.65- 




Fig.4: The total energy density for a spacial 
rotating magnetization as function of s = kR. 



Finally we want to discuss the physics of the solution and compare it with 
the experiment. 

The numerical results yield a rather short period for the modulation, 
about 3 times the radius. The reason for this short period is the smallness of 
the exchange stiffness. The sum of the demagnetization EDs {u^ + u^) /mqo 
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decreases monotonically with increasing s = kR (for constant ^o)- Only the 
exchange stiffness which increases as /c^ can hmit the value for s. 

In the real world the modulation of the magnetization has to overcome 
a serious obstacle, the pinning forces in the wire. The coercitive force is a 
manifestation of such pinning forces. In future work we intend to determine 
the strength of the nucleation force for this modulation. It has to be stronger 
than the pinning force to achieve the periodic structure. However, there are 
a number of MFM images which show a quasi-period modulation of the mag- 
netic field along the Co NW. In ref [7], Fig. 12, two MFM images are shown 
of a Co NW which is touched by a short NW. The images appear to show a 
periodic sequence of light and dark spot (in the densimeter trace along the 
NW does not resolve the fine structure). In ref. [S] the MFM image of a 
Co NW with 2R = 35nm shows a quasi-periodic field. However, the ratio of 
period to radius is not easily extracted from these images. One particularly 
good example are the experiments by Belliard et al. [2] with [Co/Cu] NWs. 
For example MFM images of a multi-wire with [170nm Co/lOnm Cu] ap- 
pear to show opposite magnetization for neighboring segments. We expect 
that the demagnetization ED causes an anti-ferromagnetic coupling between 
neighboring Co segments. 

It throws some additional light on the physics of the modulated magneti- 
zation if one applies the above considerations to a Co wire with a macroscopic 
radius, for example 2R = .8mm. If one assumes as before a modulation of 
9 = 9o cos (kz) then one obtains an optimal ED of Mt/woo = 0.213 62. This is 
about half the energy for the magnetization parallel to the wire axis. How- 
ever, for a macroscopic wire one should replace the sinusoidal phase modula- 
tion by a more favorable one, close to a square wave. This will reduce the ED 
even further. This calculation is in progress and will be published elsewhere. 

It is quite remarkable that we learn from the study of nanowires that 
the classical "ground state" of a macroscopic uniaxial wire is very different 
from what we thought it was. Of course, in the real world it will be very 
hard to prepare a macroscopic Co wire with sufficiently small concentration 
of pinning centers so that the magnetization can optimally align. Nanowires 
are much better suited for the observation of this modulation because they 
have fewer pinning centers. 

Abbreviations used: NW = nanowire, ED = energy density. 
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